Time Evolution of Spin Waves 
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A rigorous derivation of macroscopic spin- wave equations is demonstrated. We introduce a macro- 
scopic mean-field limit and derive the so-called Landau-Lifshitz equations for spin waves. We first 
discuss the ferromagnetic Heisenberg model at T = and finally extend our analysis to general spin 
hamiltonians for the same class of ferromagnetic ground states. 
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An important problem in statistical physics is to explain 
the macroscopic behavior of matter from first principles, 
in particular equations of motion of macroscopic observ- 
ables. There is by now an extensive literature on classical 
systems and we want to refer to a standard treatment by 
Spohn Q|. The general idea is that macroscopic equa- 
tions emerge from microscopic equations of motion (say 
Newton's or Schrodinger's) for the constituents of mat- 
ter by performing some kind of scaling limit. There are 
different types of scaling limits whose applicability is to 
be justified by the physical circumstances. 

In this paper we study equations of motion for the 
magnetization in certain magnetic samples. Long ago 
Landau and Lifshitz j^] derived such equations. In case of 
ferromagnets it says that the change of the magnetization 
vector, M(r, is simply given by rotation, i.e. 



dM{r,t) 
dt 



= M{r,t) X H,ft{r,t). 



(1) 



Here, the angular velocity is equal to an effective mag- 
netic field, Hcf{{r,t), which originates from the micro- 
scopic interactions of the spins. If they, say, interact only 
via ferromagnetic exchange interaction |^ with coupling 
function J, than Hcs{r,t) is given by 



7?eff(r,t) = / dr' J{r,r')Mir',t). 



(2) 



Expanding M{r',t) around (r,t) and considering an 
isotropic potential, J{r,r') = J{\r — r'\), HeS becomes 
in lowest order proportional to the laplacian AM{r,t), 
and (pi) simplifies to 



*M : (r, t) e R'' X R M{r,t) € R^. We set Bohr's magne- 
ton fio = which puts the constant on the right hand side 
equal to 1. 



d_ 

di 



M{r,t) = aM{r,t) x AM{r,t) 



(3) 



with 



a = i J drJ{\r\)\ 



Usually, (3) is called the Landau-Lifshitz equation. 

We are aware of only two other sources where 
a derivation of equation (^ is demonstrated but, never- 
theless, is partly based on non-rigorous arguments; some 
of them (e.g. magnon picture), however, are extremely 
important, influential and still a major challenge in this 
field. 

We will derive the more general equation by apply- 
ing what we call a macroscopic mean- field limit [|l|,||J|] . In 
this context we mention that the Vlasov equation for a 
quantum mean-field system has been proven by Narn- 
hofer and Sewell 0. One might think that a direct 
derivation of equation (||) is easier than of ([^) , but the for- 
mer does not seem to originate from a microscopic model 
Q. This can be compared with a recent work on the 
derivation of the Cahn-Hilliard equation |^ . We want to 
emphasize that we did not start with a mean-field model 
in the first place but we were forced, both for physical 
and mathematical reasons, to do so. We hope that this 
will become clear in our presentation. Admittedly, other 
approaches may directly lead to equation (||). 

Now we provide the set-up for the derivation of equa- 
tion (|^). Since magnetism is a purely quantum mechan- 
ical effect of interacting spin systems we clearly have to 
formulate a quantum mechanical model |l0| . For math- 
ematical convenience we consider a bounded sample, say 
the cube [0, 1]''. In order to describe the microscopic 
situation we enlarge the cube by a factor e~^, e > 0. Fi- 
nally, we let e j 0. K{e) = [0,e"i]'^ n Z'* represents the 
ferromagnet from a microscopic point of view. We use 
the general ferromagnetic Heisenberg-Kac hamiltonian, 
which is motivated by the mean-field structure of (||). 
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(4) 



to describe the microscopic exchange interaction en- 
ergy between spins in if(e). As usual, iJ' acts on the 
|iir(e)|-fold tensor product of the space C^'^^^. S is 
some fixed (but arbitrary) half-integer. Spin operators 
Sx = [Si, Si, Si) at lattice site x (acting on C^'^+^) sat- 
isfy the usual commutation relations 



[S*" , Sy"\ — i Sx 



We scale the interaction potential like J^y = 
e^'^J{ex, ey), which can be viewed as a Lebowitz-Penrose 
approximation |TT| to the true interaction; the factor e"^ 
in front of the sum provides us with an extensive energy 
function (see also another discussion below). Here, the 
function J on [0, l]^'' is assumed to be non-negative but 
not necessarily symmetric. 

The general idea relating the different levels of descrip- 
tion is to define the macroscopic magnetization profile 
M'^ (r, t) as the e-scaling limit of the expectation value of 
the time dependent microscopic spin operators S2{t), see 
(i). 

Thus the second main ingredient in our analysis, be- 
sides the hamiltonian (^, is to specify a class of (initial) 
states, denoted hereafter by g'. It is clear that (|l|) can 
only be correct at low temperature. As is well-known 
spins are parallel to each other in the (highest weight) 
ground state of the (isotropic) ferromagnetic Heisenberg 
model. In our states spins will remain almost parallel 
over short (compared to e^^) distances but vary over 
distances of order e~^, i.e. we locally describe T — 
equilibrium states as 



>o : ^ 



x£K(e) 



ol{S2)^M''{ex)}. (5) 



u}^ is a state on the matrix algebra over C^'^+^; in the 
case 5" = 1/2 the states uj^ are uniquely determined by 
MT itself. 

We rescale space by e but we do not rescale time. This 
particular space vs time scaling limit fits nicely with Lan- 
dau's idea of classical spin waves of large wave-length, 
since the microscopic spin-wave length goes as 0{e~^). 

In what follows we require M'^{r,0),cf){r) and J{r,r') 
to be continuous functions on the cube [0, 1]'*. We now 
show that at least for < |t| < C there are analytic 
measures dM'^ (r, t) such that the limits 

lime"* y (j){ex)uj'{S2it)) = [ <j){r) dM'' {r,t) (6) 



xeK{e) 



[0,1]'' 



for states G 5^ exist. Further, dM^{r,t) satisfy the 
Landau-Lifshitz equation (0-||). 

Our analysis is based on a series expansion of the time 
evolution of spin operators in finite volume: 



s2{t)^Y.-^K-ms2, 

n=0 

where K'^{H'^)S2 means the n-th commutator of H'^ with 
S2- We split this rt-th commutator into two terms, each 
of them can be obtained recurrentl}]^. 



K''{H')Sl = 



(7) 



X]j,2 y +1 '^^ans summation over all pairwise disjoint 
variables yi S K{e), i.e. yi ^ yj for all i ^ j with i,j = 
l,---,n+l, y„+i = {yn,yn+i), <3n+i = (<3n,an+i)i and 
= riLi^y.'- The polynomials P„ = P„[J^] are 
defined in the following way: 



(8) 




fc=i 



k— 1 times 



m—/c times 



denotes the matrix of the e-tensor, 6123 = 1. The rest 
term, i?^^, is given by 



K = J2K--'^mRl (9) 

k=0 

y2,...,S/fceif(c) ™=i 

with S*^™ = 11^^1=1 Sy,' and i?g = Rl = 0. Formulas 
(^-^) can be proven by induction. Since there are more 
e'^-factors than sums the operator norm of the rest term 
vanishes as e J, 0, 



lim||i?;|| = 0. 



(10) 



We first show that (for 7 = 1,2,3) the maps 

0^L^(0,t) :=lime'^ V cl,{ex)u;'{S2{t)) 

cXO ^ — ' 

xeK{t) 



^This is strictly true only for points yi such that the h- 
distance to the boundary of K{e), ||j/i — dK{e)\\i > n. But, 
for fixed n, we can always choose e to be small enough. 
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are continuous linear functional on the Banach space 
C([0, 1]'', II • lloo) (of continuous functions equiped with 
the maximum norm) for sufficiently small \t\ (de- 
pending on the interaction function J and on the 
spin S), which guarantees the existence of measures 
dM'^ (r, t) . This is accomplished by roughly estimat- 
ing e''J2xeK{e) l'/'M'^'('5'2W)l uniformly in e by ||0||i ■ 
||M||oo(l - 2t 11 J|U ||M||oo)-i (up to 0(e)) which holds 
at least for |i|-i > || J||oo ||A^||oo C"^ note that there 
are 2"n! contributing terms. This also shows analyticity 
for |i| < C, so we are allowed to interchange the e-limit 
with the Taylor expansion in time. 

By using the formulas (0-§|) for the commutator and 
the product property of the states uj'^, we can then iden- 
tify the limit 



lime'* V 



(/)(ex) 



Qn+l 



[0,1]" 



dr 4i{r) 



d 



n+l 



t=0 



M^(r, t) 



with the n-th derivative of the right hand side of (|^). 

Note, that we can neglect " surface" points (see second 
footnote) in the sum in ( ITTI ) close to dK{e), since they 
are of order e~'''~^^ and do not contribute in the limit. 

We certainly believe that everything is true for all 
times but we couldn't rigorously prove it. 

We now want to discuss some generalizations of the 
Landau-Lifshitz equations (1-3). First of all they do ap- 
pear and secondly, we can fully exploit our approach to 
these general cases. This also gives us insights in the rela- 
tion between symmetries of the microscopic hamiltonian 
and those of the spin-wave equation. 

Before changing the spin interaction we mention that 
we can replace the regular (hyper)cubic lattice structure 
Z'^ itself by other crystal structures such as polyatomic 
Bravais lattices. One could also immediately start with 
an infinitely extended sample instead of the cube (the 
test function (j) as well as the interaction J should then be 
Schwartz functions), or more rigorously by additionally 
performing a thermodynamic limit. 

We have just shown how the term M x AM descends 
from a 2-point spin interaction. In general, we might 
have contributions (in raising complexity) from an ex- 
terior magnetic field, M x B, magnetic-anisotropy en- 
ergy (due to the presence of an axis of easiest magneti- 
zation), (0,0, M'^) X M, and relativistic corrections like 
M X {M X AM), which determine the time derivative of 
M . We now explain how these terms emerge from gen- 
eral A^-point spin interactions. To this purpose we start 
with a hamiltonian 



N 



H 



N 



=iy 



^yi 



' ■ ^y^ 



(12) 



which includes multiple spin interactions up to a maxi- 
mal, but otherwise arbitrary order N. Whereas Akhiezer 
et al derive for instance M x M x AM by taking into ac- 
count higher orders in the Holstein-Primakoff expansion 
Q we directly consider general iV-point spin interactions. 

As before the tensor-valued functions {Jm)yl'.'.'.yZ 
slowly varying functions w.r.t. spatial coordinates and 
scale in magnitude as e(™-i)'^ such that the energy 
is extensive; we continue to write now {Jm)'yl'.'.'.yZ ~ 
J^-""(eyi,...,e2/™). 

Surprisingly, one can calculate the time evolution gen- 
erated by the more general hamiltonian in ( |l^ ) and ana- 
lyze the scaling limit applied to states in ^. Even more, 
the structure of the spin- wave equation is the same as in 
equation (|l|), only the effective magnetic field has to be 
generalized to 



E 

m — 1 



(13) 



(11) J,„+i(r, ri, . . . , r'J'^o^^-'^r.M-^ {r[,t) ■ ■ ■ M'^-{r'^, t). 



The formulas for the time evolution generated by Hj^ 
are analogous to (0j|) but more complicated in detail 
and will be omited here. 

We remark that the use of product-states as initial 
states for the Heisenberg model is perfectly justified by its 
local T — Q equilibrium properties, but in general cases 
we consider them just as convenient states, which seem 
(without proof) to be close to experimentally realizable 
initial states. 

All kinds of Landau-Lifshitz equations which originate 
from hamiltonian dynamics (and therefore conserve the 
magnetization) can now be derived by properly choos- 
ing the interaction potentials. In order to see how this 
works we now apply this to the interesting case of 3-point 
interactions. 



m = 



E 

x,y,z 



I T£\al3'Y qa c/3 C7 



(14) 



which after all is the motivation for our generalizations 
From the invariance under rotations w.r.t. spin indices 



{J^)xyz ^"'^'^ ■ It is fair to ask for 



we choose (J^^^^ 
spatial translation invariance of J3, such that (J|)j.j^^ = 
J(e(?/ — a:),e(z — x)). Further, it is reasonable that J 
has the properties J(j/, z) — —J{z, y) — J{—y, —z)- [Any 
odd function, /, provides an example by setting J(y, z) = 
/(|j/P — l-^l^).] After expanding the integrand in HcS at 
r we find the desired equation 



d_ 
dt 

with ( 



M{r,t) = aij M{r,t) x (M(r,i) x d.,djM{r,t)j 



/ dydz J {y, z) y^Zj 
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